Quantum disordered phase on the frustrated honeycomb lattice 
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In the present paper we study the phase diagram of the Heisenberg model on the honeycomb lattice 
with antiferromagnetic interactions up to third neighbors along the line J2 = J3 that include the 
point J2 = J3 = Ji /2, corresponding to the highly frustrated point where the classical ground state 
has macroscopic degeneracy. Using the Linear Spin- Wave, Schwinger boson technique followed by a 
mean field decoupling and exact diagonalization for small systems we find an intermediate phase with 
a spin gap and short range Neel correlations in the strong quantum limit (S = |). All techniques 
provide consistent results which allow us to predict the existence of a quantum disordered phase, 
which may have been observed in recent high-field ESR measurements in manganites. 
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The two-dimensional Heisenberg model in frustrated 
geometries has received a lot of attention in the last years 
One of the reasons for this interest is the common 
belief that geometrical frustration in two-dimensional 
(2D) antiferromagnets may enhance the effect of quan- 
tum spin fluctuations and hence suppress magnetic order 
giving rise to a spin liquid [B| . One candidate to test these 
ideas is the honeycomb lattice, which is bipartite and has 
a classical Neel ground state, but due to the small coor- 
dination number (z = 3), quantum fluctuations could be 
expected to be stronger than those in the square lattice 
and may destroy the antiferromagnetic long-range order 
(LRO). 

The study of frustrated quantum magnets on the hon- 
eycomb lattice has also experimental motivations: on the 
one hand, recent ESR experiments in high magnetic field 
on Bi 3 Mn 4 0i2(N0 3 ), which is described by a S = 3/2 
honeycomb Heisenberg model, have lead to the conjec- 
ture that geometric frustration plays an important role 
in removing the long range magnetic order Q. On the 
other hand, the family of compounds BaM2{XO/±)2 with 
M = Co,Ni and X = P, As, consists of magnetics ions 
M arranged in weakly coupled frustrated honeycomb lat- 
tices with spin S = 1/2 for Co and S = 1 for Ni 0. 

Last but not least, new possibilities to design inter- 
acting lattice systems with controlled geometry emerge: 
modern strategies in chemistry open a route to synthe- 
size new materials with a desired lattice structure and 
intersite interaction Q and the controlled setup of op- 
tical lattices for cold atoms would allow to create arbi- 
trary lattice structures as well as to tune the interactions 

In the present paper we study the Heisenberg model 
on the honeycomb lattice with first (Ji), second ( J2) and 
third (J3) neighbors couplings [HJ, along the special line 
<h = Jz- Using Linear Spin- Wave theory, the Schwinger 
boson technique and exact diagonalization we find strong 
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FIG. 1: (Color online) Classical phase diagram for the model 
([T]). There are three different phases schematized by the car- 
toons: The Neel phase (I) has two (antiparallel) ferromagnetic 
triangular sublattices (blue and red balls) , the collinear phase 
(II) where ferromagnetic chains are antiparallel and the spi- 
ral phase (III). We focus on the dashed line corresponding to 
J2 = J3 that includes the special point (J2 = J3 = Ji/2) 
where the ground state is infinitely degenerate. 



evidence for the existence of an intermediate disordered 
region where a spin gap opens and spin-spin correlations 
decay exponentially. Although our results correspond to 
a specific line, we conjecture that the quantum disordered 
phase that we have found in the vicinity of the tricritical 
point extends in a finite region around it. 



The Heisenberg model on the honeycomb lattice is de- 
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scribed by 

H = Ji ^2 Si ■ S 3 ,+ J 2 ^2 Si • + J z J2 



were, Si is the spin operator on site i 
The classical model 
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12j displays different zero tem- 
perature phases shown in Fig. (QJ with a tricritical point 
at J2 = J3 = At this point, the Hamiltonian can 

be written, up to a constant term, as a sum over edge- 
sharing plaqucttcs 
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where S 7 = X^e 7 ^ is the total spin in the hexagon 7. 

At this particular point any state with vanishing total 
spin per elementary hexagon is a classical ground state, 
giving rise to a large GS degeneracy. This kind of sit- 
uation is reminiscent of what happens in the J\ — J% 
square lattice case for J 2 /,h = 0.5 [H, 

which presents 

a disordered spin gapped phase around this point at the 
quantum level. 

Motivated by this analogy, we study here the region 
around the tricritical point Ji = J3 = \J\ keeping Ji = 
J3 as shown in Fig. [1] 

We construct in the following the corresponding mag- 
netic phase diagram with the results obtained by using 
spin-waves, Schwinger bosons and ED techniques. 

The classical phase diagram reduces to that shown in 
the line 1/S C = of Fig [2] where two collinear phases 
meet at the classical critical point Jij J\ = 0.5. A general 
ordered planar spin configuration can be written as: 



4a =S(C0S [Q-R 



ei 



sin 



(g-i? + Q )e 2 )(3) 



where R is the position of the unit cell and ei , ei are the 
primitive vectors of the direct lattice. The magnetic con- 
figuration is characterized by the ordering wavevector Q 
and the internal phase <p a , where a = A, B is an internal 
index on each sublattice. There are two different phases 
for the classical ground state: for Jij J\ < 0.5 the ground 
state corresponds to the Neel phase, with ordering vec- 
tor Q = (0,0) and 4>a — 4>b = n and for Jij J\ > 0.5 
the ground state corresponds to a collinear phase with 
Q = (27r/\/3, 0) and 4>a — 4>b = tt. In the following we 
call this phase columnar antifcrromagnetic (CAF) phase. 
Using linear spin-wave theory (LSWT) HE!, we have 
investigated the stability of the classical phase diagram. 
The results obtained are summarized in Fig. [5] as a func- 
tion of the spin-S 1 and the frustration (J1/J1) (we do 
not present details in the different steps of LSWT since 
these are standard 14, IE 12 Q- The edge of the ordered 
phases was obtained by finding the frustration ( J2/J1) at 
which the quantum fluctuations destroy the classical or- 
der, i.e. where the order parameter vanishes. The results 
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FIG. 2: (Color online) Inverse of the critical pin S c as a 
function of J2: The dashed-point line shows the bound ob- 
tained by LSWT. Above this line, quantum fluctuations de- 
stroy the magnetic order, the staggered magnetization go to 
zero. For 5=1/2 (full black line), there is a small range 
0.29 < J2/J1 < 0.55 where there is not magnetic order; the 
full blue-line corresponds to the edge obtained by SBMFT. 
For the 5=1/2 case, there is a range 0.41 < J2/J1 < 0.6 
where the system has a spin-gap indicating a quantum dis- 
ordered phase (see Fig. [4j. The dotted- line correspond to 
the classical limit S — > 00 where the ground state corre- 
spond to the Neel phase in the region J2/J1 < 0.5 , while for 
J2/J1 > 0.5 the ground state correspond to the CAF phase. 



show that quantum fluctuations reduce the Neel stabil- 
ity around the classical point Jij J\ = 1/2. For S = 1/2, 
Neel order is found for < J 2 /Ji < 0.29 and the CAF 
order for 0.55 < Ji/J\- To further analyze the region 
0.29 < Ji/ J\ < 0.55 which is seen as a non-magnetic 
phase within LSWT we now use the Schwingcr-boson 
mean field theory (SBMFT) and Exact Diagonalization. 

The Heiscnberg model on the honeycomb lattice was 
studied using SBMFT by Matson et al [11] for antiferro- 
magnctic interactions at first and second neighbors. Here 
we study Hamiltonian ([1]) using a rotationally invariant 
version of this technique, which has proven successful in 
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incorporating quantum fluctuations pi 

In the Schwinger-boson approach, the Heisenberg in- 
teraction is written as a biquadratic form. In this repre- 
sentation the spin operators on each site are replaced by 
two species of bosons via the relation: 



Sg = -ht.a.hg, 



(4) 



where bt = (b]L, bt.) is a bosonic spinor, a is the vector 
of Pauli matrices, and there is a boson-number restriction 
Sct bt^bgo- = 25 on each site. 

With this representation, the rotational invariant spin- 
spin interaction can be written as 



(5) 



where a = zb, : O : indicates the normal ordering of the 



3 



operator O and we define the SU (2) invariants 

A .'.,/ = ^^^br.„b ;/ . „ (6) 
Bt- = -Vbt b 9a . (7) 

cr 

The operator Af.j creates a spin singlet pair between 
sites x and y and ~Bs,y creates a ferromagnetic bond, 

I 



which implies the intersite coherent hopping of the 
Schwingcr bosons. This rotational invariant decoupling, 
which enables to treat ferromagnetism and antifcrromag- 
netism on equal footing, was introduced by Ceccato et al 
[l9j| . Later, Flint and Coleman [2(| presented a general- 
ization to large N that shows substantial improvements 
over the Sp(N) approach. 

To generate a mean field theory, we perform the 
Hartree-Fock decoupling 
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[Bapi®- y) B S+r«,y+r 3 ~ A* a g (x - y) A S+ ? a £ +Ff> + H.c] - ((Saf+f^ ■ S^+p ? )mf) 



(8) 



where 

A* aS (x-ff) = (Al +r ^ +P/3 ) (9) 

((Sx+r a ■ Sy+r l3 )MF) = \B aj3 {x - y)\ 2 - | A aB (x - y ) 1 2 

(10) 

These are the mean field equations and must be solved 
in a self-consistent way together with the constraints on 
the number of bosons 

B aa (R = Q) = 2N C S. (11) 

where N c is the total number of unit cells and 5* is the 
spin strength. Numerical solutions of the ([9]) and (fTTj) 
involves finding the roots of 24 coupled nonlinear equa- 
tions for the parameters A and B plus the additional 
constraints to determine the values of the Lagrange mul- 
tipliers X^K We perform the calculations for finite but 
very large lattices and finally we extrapolate the results 
to the thermodynamic limit [24) . 

We solve numerically ([9]) and (|lip for several values 
of the frustration parameter J2/J1 and with the values 
obtained for the MF parameters and the Lagrange mul- 
tipliers we compute the energy and the new values for 
the MF parameters. We repeat this procedure until the 
energy and the MF parameters converge. After reach- 
ing convergence we can compute all physical quantities 
like energy, spin-spin correlations and the excitation gap. 
In order to support the analytical results of the MF ap- 
proach, we have also performed exact diagonalization on 
finite systems with 18, 24 and 32 spins with periodic 
boundary conditions for 5 = 1/2 using Spinpack [23| . In 
Fig. EJa) we show the ground state energy per unit cell 
as a function of the frustration for a system of 32 sites 
calculated by mean a SBMFT and ED showing an excel- 
lent agreement between both approaches. The advantage 
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FIG. 3: (a): GS energy per unit cell as a function of J2/J1 
for a lattice of 32 sites. The circles are exact results (ED) and 
the squares are the SBMFT results, (b): Spin-Spin correla- 
tion function (SSCF) vs distance X in the zig-zag direction 
obtained within SBMFT. For < J2/J1 < 0.41, the SSCF 
correspond to the Neel phase white long-rage-order (LRO), 
for 0.41 < J2/J1 < 0.6 the correlations are short ranged indi- 
cating a gap zone with sort-range-order (SRO). 

of the SBMFT is that it allows to study much larger sys- 
tems: we have studied different system sizes up to 3200 
sites, which we extrapolate to the thermodynamic limit. 

For the present model we only find commensurate 
collinear phases (CCP) and for these phases the wave 
vector Qo, where the dispersion relation has a minimum, 
remains pinned at a commensurate point in the Bril- 
louin zone, independently of the value of the frustration 
J2/J1. In the thermodynamic limit N c — > 00, a state 
with LRO is characterized in the Schwinger boson ap- 
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FIG. 4: (a): Gap in the boson dispersion as a function of 
J2/J1 for S — 1/2. In the region J2/J1 ~ 0.6 the system 
remains gapped. Inset: finite size scaling for the gap. (i) 
J2/J1 = 0.5 (7 = 0.6451), (ii): Circles correspond to J2/J1 = 
0.05 (7 = 0.911) and squares correspond to J2/J1 = 0.35 
(7 = 0.758). 



for 5 = 1/2. On the one hand, LSW analysis predicts 
that the Neel LRO region extends up to J2/J1 ~ 0.29 
and the collincar antifcrromagnct LRO is present for 
J2/J1 ^ 0.55. In the intermediate region we find no evi- 
dence of any type of magnetic order with this technique. 
On the other hand the results found with SBMFT pre- 
dict a quantum disordered region 0.41 < J2/J1 < 0.6. In 
this region a gap opens in the bosonic dispersion and the 
analysis of the spin-spin correlation function shows Neel 
short range order followed by the LRO CAF phase for 
J2/J1 > 0.6. 
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proach by a condensation of bosons at the wave vector 
Qq. This implies that the dispersion of the bosons in 
a state with LRO is gapless. As we discussed earlier, 
we solve ([9]) and (jTTJ) for finite systems, then to detect 
LRO we calculate the gap in the bosonic spectrum as 
a function of J2/J1 for different system sizes and per- 
form a finite size scaling finding a finite region where 
the system remains gapped. The extrapolation of the 
boson gap as a function of the frustration is shown in 
Fig. |U In the range of 0.41 < J2/J1 < 0.6 the sys- 
tem presents a finite spinon gap. The inset shows an 
example of the finite size scaling for different values of 
the frustration. The structure of the different phases can 
be understood calculating the spin-spin correlation func- 
tion (SSCF). For J2/J1 < 0.41 the SSCF is antiferro- 
magnetic in all directions and is long-ranged while for 
0.6 < J2/J1 we have found ferromagnetic LRO corre- 
lations in the zig-zag direction that correspond to the 
CAF phase. The most interesting result is in the inter- 
mediate region 0.41 < J2/J1 < 0.6 where the results for 
the SSCF are consistent with the presence of a gap. In 
this region we have found short-ranged antiferromagnetic 
correlations. A plot of the SSCF for J2/J1 = 0,0.2 and 
0.6 obtained within SBMFT is presented in Fig. [3b. 

The Fig. [2] shows the ground state phase diagram as 
a funtion of 1/5. On one hand, for 1/5 smaller than a 
critical value l/5c(J2/Ji), the correlation function have 
LRO, characterized by a condensation of bosons at the 
wave vector Q Q . On the other hand, when 1/5 greater 
than l/5c(J2/Ji), the correlation functions have SRO 
indicating quantum disorder. 

In summary, the results obtained with the different 
techniques used in the present paper suggest the exis- 
tence of a region in the intermediate frustration regime 
where the system does not show quantum magnetic order 
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